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Abstract. The Fourier transforms of Laguerre functions play the same canon- 
ical role in wavelet analysis as do the Hermite functions in Gabor analysis. We 
will use them as analyzing wavelets in a similar way the Hermite functions were 
recently by Grochenig and Lyubarskii in Gabor frames with Hermite functions, 
C. R. Acad. Sci. Paris, Ser. I 344 157-162 (2007). Building on Seip's work 
Beurling type density theorems in the unit disc, Invent. Math., 113, 21-39 
(1993), concerning sampling sequences on weighted Bergman spaces, we find 
a sufficient density condition for constructing frames by translations and dila- 
tions of the Fourier transform of the nth Laguerre function. As in Grochenig- 
Lyubarskii theorem, the density increases with n, and in the special case of 
the hyperbolic lattice in the upper half plane it is given by 

4tt 

ologa < , 

2n + a 

where a is the parameter of the Laguerre function. 



1. Introduction 

The "frame density problem" is one of the fundamental questions in applied 
harmonic analysis: how coarse we need to sample a continuous object so that the 
resulting discrete object is a frame in a certain Hilbert space? This question was 
sharply solved for frames in the Bargmann-Fock [8] , [24] , |25j and in the Bergmann 
space [28 . However, in the cases of general wavelet and Gabor transforms, very 
little is known and only a few very special windows and analysing wavelets are 
understood. 

The short time Fourier (Gabor) transform with respect to a gaussian window 
can be written in terms of the Bargmann transform, mapping isometrically the 
space L 2 (R) onto the Bargmann-Fock space of entire functions. This is the reason 
why everything is known about the geometry of sequences that generate frames by 
sampling the Gabor transform with gaussian windows g(t) — e _7r * . Apart from 
this example, the only cases where a description is known of the lattice sequences 
that generate frames are the hyperbolic secant g(t) = (cosh at) ~ 1 [18j and the 
characteristic function of an interval [19] , which turned out to be a nontrivial 
problem. There is also a necessary density condition for Gabor frames due to 
Ramanathan and Steger |23] and it was recently observed that the frame property 
is stable under small perturbations of the lattice, given the window belongs to the 
Feichtinger algebra [10] , 
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If we consider the wavelet case, things appear to be even more mysterious. Al- 
though the density and irregular grid problems have attracted some attention lately 
(see pQ, [31] . [22] . |15j and references there in), there is no known counterpart of 
Ramanathan and Steeger theorem (there is, however, a version of the HAP property 
for wavelet systems |16j ) and the only information available so far, on particular 
systems, concerns a special family of analyzing wavelets that maps the problem into 
Bergmann spaces: The wavelet transform, with positive dilation parameter, with 
respect to a wavelet of the form (Paul's Wavelet in some literature) 

can be rescaled as an isometrical integral transform between spaces of analytical 
functions, namely, between the Hardy space on the upper half plane H 2 (TJ) and 
the Bergman space in the upper half plane ^4 2q+1 (U). It is natural to refer to 
this isomorphism as the Bergman transform (the designation analytic wavelet is 
also frequently used). Since the upper half plane can be mapped isomorphically 
onto the unit disc by using linear fractional transformations, we can construct a 
transform mapping the Hardy space onto the Bergman space in the unit disc. 

The approaches used to deal with the special situations mentioned in the above 
paragraphs are based in techniques which differ from case to case. It seems highly 
desirable to follow a more structured approach. The natural place to look for this 
structure is within the context of Hilbert spaces of analytic functions, where the 
powerful methods from complex analysis may answer questions that seem hopeless 
otherwise. Following this line of reasoning implies using windows that allow to carry 
the problem to such spaces. In this direction a major step was taken recently by 
Gochenig and Lyubarskii [12j . by considering Gabor systems with Hermite func- 
tions of order n as windows. They have proved that if the size of the lattice A is 
< (n + 1) _1 then the referred Gabor system is a frame and provided an example 
supporting their conjecture that the result is sharp. 

In the Bargmann-Fock setting the Hermite functions play a very special role 
[111 pag. 57]. They are, up to normalization constants, pre-images, under the 
Bargmann transform, of the monomials {z n } and since the latter constitute an 
orthogonal basis of the Bargmann-Fock space, their pre-images also constitute an 
orthogonal basis of L 2 (R). That is, to the isomorphism 

(1.2) L 2 (R)^F 2 (C) 
corresponds 

(1.3) h n ^ c n z n 

where B is the Bargmann transform, F 2 (C) is the Bargmann-Fock space, h n are the 
Hermite functions and c n some constants dependent on n. They are canonical to 
time-frequency analysis in an additional sense, since they constitute the eigenfunc- 
tions of the time- frequency-localization operator with Gaussian window [3]. The 
Hermite functions are eigenfunctions of the Fourier transform and they can be used 
[20] to describe Feichtinger 's algebra S . 

Wavelet and Gabor analysis share many similarities and many of their struc- 
tural aspects can be bound together in a more general theory using representations 
of locally compact abelian groups [13] . [9]. In looking for a Wavelet-analogue of 
Grochenig-Lyubarskii structured approach to the density problem, we must first 
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clarify what functions should be used instead of the Hermite functions. In analogy 
to the above paragraph it is natural to consider the pre-images, under the Bergman 
transform, of the monomials {z n } in the unit disc (up to an isomorphism with the 
upper half plane). Since the {z n } form a basis of the weighted Bergman spaces in 
the unit disc (independently of the weight), their pre-images must be an orthogonal 
basis of the Hardy space H 2 (U) and it is reasonable to expect that such functions 
will play a similar role in Wavelet analysis as do the Hermite functions in Gabor 
analysis. Such functions are the Fourier transforms of the Laguerre functions and 
they will be ad-hoc denoted by S" . 

Additional structural evidence that the functions S" are canonical in Wavelet 
analysis comes from the work of Daubechies and Paul [4] , where it is shown that 
they are the eigenf unctions of a differential operator that commutes with a time- 
scale localization operator, once windows of the form (|1.1|> are chosen. This com- 
pletely parallels the situation with the time-frequency-localization operators with 
gaussian windows, which commute with the harmonic oscillator and therefore have 
as eigenfunctions the Hermite functions [3] . It was observed by Seip |26j that these 
problems have a more natural formulation when mapped into the convenient spaces 
of analytical functions. 

2. Description of the results 

There exists an analogue of the correspondence between (jl.2p and (|1.3[) . but 
involving four different functional spaces and the corresponding bases: To the se- 
quence of isomorphisms between the Hilbert spaces 

(2.1) L 2 (0,oc) H 2 (U) B ^ A 2a+1 (V) ^ A 2a+1 (T>), 

where A2 a +i(U) and ^2a+i(D) are the weighted Bergmann spaces in the upper- 
half plane and in the unit disc, respectively, corresponds the relations between the 
basis of the respective spaces: 

(2.2) c s: B z a cX ^ <z n . 

for some constants c", where I™ is the nth Laguerre function of order n and pa- 
rameter a, S" is its Fourier transform and 9" is a basis of A 2a +i(U) to be defined 
in section 4. This correspondence was implicit in the conection between papers [4] 
and [26] but since it was not stated explicitly, we devote section 4 to clarify how 
exactly does it work. 

The functions S" were computed recently in closed form [29] , but the connec- 
tion to the Wavelet transform seems to have been unnoticed. It is also shown in 
[29] that S" are, up to a fractional transformation, defined in terms of a certain 
system of orthogonal polynomials on the unit circle. We will show that this family 
of polynomials is nothing more but the circular Jacobi orthogonal polynomials for 
which there is, for computational purposes, a very convenient three term recurrence 
formula. This connection will make our case to call the functions S" the rational 
Jacobi orthogonal functions. 

With a computational method available for evaluating the functions S" and 
graphic evidence of their good localization properties (see for example the plots of 
their real versions in page 44 of [2]), it is natural to investigate how to use such func- 
tions as analysing wavelets and to obtain frames from the resulting discretization. 



4 



LUfS DANIEL ABREU 



We will obtain the following sufficient condition on the density of the parameters 
of the hyperbolic lattice {(a^bk, a J )}j t kez' 

&loga < - — ■ — 
2n + a 

that is, as in Grochenig and Lyubarskii result [12j . this density increases with n. 
This will follow as a special case of the following more general theorem, which uses 
Seip's notion of lower Bcurling density in the unit disk and constitutes the main 
result in this paper: 

Theorem 5.1 Let Y C D denote a separated sequence obtained from mapping 
the sequence {z m j — a m ,j + ib m ,j} C U into the unit disk via a Cayley transform. 
If D-{T) > n + f then { !„., Ih„ S£(f)} mJ is a frame of H 2 (U). 

Case a = gives the sufficient part of theorem 7.1 in [28] . The proof of theorem 
1 will be dramatically simplified with the observation of the following remarkable 
ocurrence: There exist constants {C} and sequences of numbers {a^}, both de- 
pending on a and n, such that 

t - 

(2-3) S^(-)=Cj2^ k+q (t). 

This will allow to obtain a formula which expresses the wavelet transform with 
window S" as a linear combination of simple but non-analytic functions. As a 
result we will have to deal with a situation that is reminiscent of the non-analiticity 
problem that Grochenig and Lyubarskii faced with formula (15) of [12] and which 
they solved by using the so called Wexler-Raz identities |llj . Since there is no 
known analogue of the Wexler-Raz identities in wavelet analysis, we have to deal 
with this in a more direct way. The key idea is to map the problem into the unit 
circle and use a deep sampling formula proved by Seip in [28] , Althought this 
approach lacks the simplicity of [12] it gives the result for more general sequences, 
while the method in [12] works only for the lattice case due to the restritions 
imposed by the use of the Wexler-Raz identities. 

We organize our ideas as follows. The next section contains the main definitions 
and facts concerning wavelet transforms, Bergman spaces and Laguerre functions. 
Section four contains the evaluation of the pre-images required to build up the 
correspondences (|2.1j) and (|2.2|) . The fifth section contains our main results on 
wavelet frames with Fourier transforms of Laguerre functions. We conclude the 
paper collecting some further properties of the functions S" and clarifying their 
classification within known families of orthogonal polynomials. 

3. Tools 

3.1. The Bergman transform. Now we present a sinthesis of ideas that 
appeared in the section 3.2 of [13] and in [4] (see also [2j pag. 31]) Here they will 
be exposed in such a way that the role of the Bergman spaces is emphasized. 
Consider the dilation and translation operators 

D s f(x) = isrhis- 1 *) 

T x f(t) = f(t-x) 

and define 

r' s {t) = T x D s ^{t) = \s\-^ ). 
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The wavelet transform of a function /, with respect to the wavelet tp is 

/oo 
-oo 

A function ip 6 L 2 (R) is said to be admissible if 

/ \T^{s)\ 2 - = K 
Jo s 

where if is a constant. If ip is admissible, then for all / <E L 2 (R) we have 

poo poo 

(3.1) / / s- 2 \W^f(x,s)\ 2 dxds = K\\f\\ 2 



We will restrict ourselves to parameters s > and functions / 6 H 2 (U) , where 
H 2 (XJ) is the Hardy space in the upper half plane U = {z = x + is : s > 0} 

/oo 
\f(x + is)\ 2 dx < oo}. 
-oo 

Let T denote the Fourier transform 

^ poo 



(Ff)(t) = -j= I e-^f(x)dx 

V ZTt J-oo 

By the Paley- Wiener theorem, H 2 (U) is constituted by the functions whose Fourier 
transform is supported in (0, +oo) and belongs to L 2 (0,oo). 
Now take the special window ip a (t) defined in Since 



(3.2) = l [0 , oo] s a+ h a e i ^ + ^ t 

then 



(3.3) W^J(-x,s) = s a+ 2 / t a e lzt {Tf){t)dt 

Jo 

where the function defined by the integral is analytic in z = x + si. The identity 
(|3.1|) gives 



(3.4) 



J J \ W i> a f(- x ' s )\ 2 s 2 dxds = ll/H ff2(R) 



This motivates the definition of the Bergman transform, or the analytic wavelet 
transform: 



(3.5) Ber c 



p OO 

f(z)= / t a e izt (Ff)(t)dt = s- a -iWT P J(-x,s), 
Jo 



where z = x + is (see for instance [15] . where the authors use this Bergman trans- 
form with the same normalization in the case a = 1). Introducing the scale of 
weighted Bergman spaces 



A a (TJ) = {/ analytic in U such that / / |/(z)| 2 s Q 2 dxds < oo}, 



it is clear from (|3 .4[) and (|3.5[) that Ber a f(z) € A2 Q +i(U). We have therefore an 
isometric transformation 

Ber a :iJ 2 (R)^ A 2a+1 (U) 
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The weighted Bergman spaces in the unit disc are denoted by ^4 a (D) and defined 

as 



A a (D) = < f analytic in D such that 



\f(z)\'(l-\z\) a -'dxdy<oo 



For detailed expositions of the theory of Bergman spaces we point out the mono- 
graphs m, M, M 

3.2. Frames and sampling sequences. A sequence of functions {ej} is said 
to be a frame in a Hilbert space H if there exist constants A and B such that 



(3.6) 



A\\f\\ 2 <^|(/, e ,)| 2 <i3 ll/H 



In spaces of analytic functions a related concept to frames is the one of a 
sampling sequence. A set T — {zj} is said to be a sampling sequence for the 
Bergman space A a (U) if there exist positive constants A and B such that 



(3.7) A 



u 



\f(z)\ 2 y a - 2 dxds<Y,\f(zj)\ 2 y?<B 



\f(z)\ 2 y a - 2 dxds 



u 



The corresponding definition of a sampling sequence for the Bergman space on 
the unit disc is 

a[ [ \f(z)\ 2 (i-\z\ 2 r-*dz < ^2\f( Zj )\ 2 (i-\z\ 2 r 



< B 



\f(z)\ 2 (i-\ z \ 2 r- 2 dz. 



Define the pseudohyperbolic metric on the unit disk by 



l-(z 



> and its 



Following [28 , a sequence T — {zj} C D is separated if inf ra ^j 
lower density D~ is given by 

D"(T) = lim inf inf 1 ^ " )] 

r^l z log ^ 

The next Theorem is from [28j : 

Theorem A A separated sequence T C D is a sampling sequence for A a (D) 
iff D-(T) > | - i. 

It is a generalization of an earlier lattice result |27j : 

Theorem B Let T(a,b) = {zjk}j t kez, where zj^ — a J (bk + i). T(a 1 b) is a 
sampling sequence for A a (TJ) if and only if &loga < -^zj- 

3.3. Fourier transforms of Laguerre functions. The Laguerre polynomi- 
als will play a central role in our discussion. One way to define them is by means 
of the Rodrigues formula 



(3.8) 



n\ dx r 



[e- x x a+n ] 
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and this gives, in power series form 

T o.( _ ( a + !)» i- n )k x k 

For information on specific systems of orthogonal polynomials, we suggest |17j . The 
Laguerre functions are defined as 

l«(x) = l [0 , oo] ( X )e-^ 2 x^L«(x) 

and they are known to constitute an orthogonal basis for the space £ 2 (0, oo). By the 
Paley- Wiener theorem, the Fourier transform is an isomorphism between H 2 (JJ) 
and L 2 (0,oo). Therefore, the Fourier transform of the functions I" form an orthog- 
onal basis for the space H 2 (XJ). From now on we will set 

(3.9) TS«{t) = l°(t) 

The functions S"(t) can be evaluated explictly using (|3.8[) . This was already done 
by Shen in |29j , where the author was interested in describing the orthogonal poly- 
nomials which arise from an application of the Fourier transform on the Laguerre 
polynomials. A description of this method of generating new families of orthogonal 
polynomials with an interesting historical account is in the paper |21j where a link 
is provided between Jacobi and Meixner-Polaczek polynomials. Following [29], we 
have 



hi — ' 



k=Q k\(a+l) k \±-it 
4. Bergman transform of S"(t) 



In this section the correspondences (|2.1j) and (|2.2p will be established via a 
direct calculation. The complex linear fractional transformations will play an im- 
portant role, in a style that is reminiscent of the way they are used in the discrete 
series representation of SL(2, R) over the Bergman space chapter IX]. 

We first define a set of functions that, as we shall see later, constitute a basis 
of A a (XJ). For every n > and a > —1 let 

1 \ n / in -a— 1 

2 > ' ■ 1 



IZ 2 



Now define a map T a such that for every function / 6 A a (U) the action of T a 

is 

a + l 

j 1 1 — ]— i / i \ 

(4.1) T a f(w) = f 



2w 



1 — w 



The range space of T a is a weighted Bergman space in the unit disc. 
Lemma 4.1 The map 

T a . A a (TJ) ► A a (D) 

is an unitary isometry between Hilbert spaces. 

Proof. Argue as in the proof of lemma 1 in [3 pag.185]. □ 
Proposition 4.1 For n = 0, 1, the following relations hold: 

(4.2) Ber a {S 2 n a ) = « , with < = (-l) tt+1 (2a + n )!/n! 
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and, for \z\ < 1, 

(4.3) = z n 

In other words, the function S 2a is the pre-image, under T a o ^-^Ber a ^j , of z n € 
A a (D). 

Proof. Using the definition of the Bergman transform (|3 . 5[) and (|3.9p we have 

/>oo 

Ber a S 2 n a {z) = / t a e lzt (FS 2a (t))(t)dt 







e (i*-h)t t 2a L 2 n a (t)dt 



dx n 



d n 



where in the last two identities we are applying Rodrigues formula (|3.8[) for the 
Laguerre polynomials and writing the integral in terms of the Laplace transform C, 
whose well known properties establish the last identity. Now, the linear fractional 
transformation 



is an analytic isomorphism between the upper half plane and the unit circle. Since 
the inverse of this transformation is given by 



i w + 1 



z = — - 



2w-r 

a short calculation with the definition of T a gives (|4. 31) . □ 

Remark 1. {^"(z)} is a basis of A a (U) and the map T a is an unitary iso- 
morphism between A a (U) and A a (D). Indeed, since {z n } is an orthogonal basis of 
the space A a (D) pag. 186] and it is contained in the range of T a , T a is onto 
and therefore an unitary isomorphism. The functions {^"(z)} form a basis for the 
space A a (\J) since they are the pre-images of the basis {z n }. 

As a consequence we obtain a new proof of the (known) isomorphic property 
of the Bergman transform. 

Corollary 1. The transform Ber a : H 2 (R) — » A2 Q +i(U) is an isometric 
isomorphism. 

Proof. The isometry is a consequence of the isometric property of the wavelet 
transform, so we need only to prove that the Ber% is onto. But in view of the 
preceding remark and Theorem 5.1, the range of Ber% contains a basis of A a (TJ). 
Therefore, Ber^ is onto. □ 

Remark 2. Similar calculations as we have seen here also play a role in [5], 
in the context of Laplace transformations and group representations and in [4], to 
obtain an explicit formula for the eigenvalues of the time-scale localization operator. 
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5. Wavelet frames with Fourier transforms of Laguerre functions 

5.1. Preparation. We wish to construct wavelet frames with analysing wavelets 
5™, by using the common discretization of the continuous wavelet transform via 
the hyperbolic lattice. Discretizing the scale parameter s in the Wavelet transform 
by a sequence a J and the parameter x by a J bm gives 

(5.1) W^f{o?bm, a 1 ) = (/, T a3bm D a ^) . 

We want to know conditions in a and b under which {T a j^D a jip} is a frame, for a 
given wavelet ip. More generally we can consider the non lattice problem replacing 
(a J 6m, a J ) by more sequences (aj m ,bj m ). 

Remark 3. Observe that it follows immeadiately from Theorem A and Theorem 
B that ifT C D is the sequence obtained from mapping the sequence {Zk,j — a J (bk + 
i)} C U into the unit disk via a Cayley transform, then D~(T) = b j^ a (to move 
between the notation of [28.J and [2? set a = In + I). 

Remark 4. Observe that {T a j bm D a jip a }j t k^z is a frame ofH 2 (XJ) iffT(a,b) = 
{z = a J bm + a J i}j^kez is a set of sampling for A2q.+i(U). Indeed, since functions 
in A2a+i(U) can be identified with Ber a transforms of H 2 (XJ) functions, it follows 
from \3. 7| ) and (?? ) that T(a, b) is a set of sampling for A2 Q +i(U) if and only if 

A[[ \W^J(x, S )\ 2 ^<J2\W^J(a^b m ,a^\ 2 <B J J \W+J{x, stf 

U j,m U 

Using \3.4-\ l this is equivalent to 

A ll/H H 2 (V) 

which says that {T a j bm D a jip a }j t k is a frame in iJ 2 (U). 

We now collect the required preliminary lemmas for the proof of the main 
result. The first was quoted in the introduction and is just a simple modification 
of the representation (|3.10[) . The second appeared in [27] (see [6j pag. 160] for a 
proof) and the third one is a quite deep sampling theorem which is formula (30) in 
[28] with s = and e = (see also [6j pag. 216] for a direct derivation with this 
choice of parameters). 

Lemma 5.2 The functions 5 1 " can written as the linear combination 12. 3]) of 
analyzing wavelets ip k+ a(t) defined by 

t " 
S? t (-)=CY,^ k+§ (t), 

where the constants C and the coeficients a k given as 

T(f + 1) (l + a) n 



!(-n) fc (f + l) fe 



k\{a + l) k 
Lemma 5.3 For 1 < s < t we have 

(5.2) E ^'-'T ^ga-ici 2 )'-* 

.7 1-C«J 
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for any separated sequence {zj}. 

Lemma 5.4 Let D({wj}) > (3. Then every f 6 A2p+i(D) satisfies 

(5.3) (1 - MY/H - £(1 - K-|Y/K)ft», 

3 

with the estimate for hj (w) 

(5 . 4) IM.)|< C "-H'>ii-'"A 

|1 - WjW\ 

5.2. Main result. To prove our main result we will adhere to the following 
plan: we begin by deriving from Lemma 5.2 a formula expressing the wavelet trans- 
form with analysing wavelet in terms of the wavelet transform with analysing 
wavelet ip k+ a. The proof of the right frame inequality is then straightforward (as 
usual in frame inequalities). To prove the left hand inequality we first rewrite the 
formula in terms of complex functions defined on the upper half plane and then map 
it into the unit disk by means of a Cayley transform. Manipulating the resulting 
formula and using lemma 5.3 and lemma 5.4 gives an estimate in the unit disk. We 
proceed backwards to the half plane to end up with the required estimate in terms 
of wavelet frames. 



Theorem 5.1 Let V C D denote a separated sequence obtained from mapping 
the sequence {z m j = a mj + ib m j} C U into the unit disk via a Cayley transform. 
Lf D-(T) > n+ § then '{T amJ D bnhi S%(±)} md is a frame of H 2 {U). 

Proof. The definition of the wavelet transform gives, using (|2.3p . 



,T_ X D S S%(-) 



W ss{i) f(-x,s) = (/,'. 

n 



k=0 



C^OkWv, k+f f(-x,s). 



k=0 



The right frame inequality is now easy: observe that, since {T am j Db m jip^+f} 
is a frame, then for every k = 1, n there exists a B k such that 

2 



E <B k \\f\\ H 2 (U) ■ 



Let B n = maxKAKn B k . Then, 



E 



f,T amd D bm ^{-) 



E 



k=0 



< Cj2\ak\ 2 J2\(f> T ^ D ^ 



fe=0 



< B\\f\\ 2 , 
where B = B n C^Jl =0 \ak\ 2 and we are done. 
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Now the left frame inequality. Define the function F(z) as 

n 

(5.5) F(z) = W Ss(i) f(-x, s) = Cj2akS k+ % +L *Ber k +%f(2 



k=0 



where z = x + is belongs to the upper half plane U. We construct a related function 
in the unit disc, by setting w — jt^' 

n 



k=0 



where 



and 



1 — w 



Gfc+f (w) 



1 — w 



2k+a + l 



,1 + W 



Ber k+ *f{i-—) G ,4 2fc+Q+1 (D). 
1 — w 

Therefore, if D({wj}) > n + ^ we can use (|5.3p to write, for every k < ) 
(1 - \w\ 2 ) k+ %G k+ «(w) = £(1 - Kf) fc +t G fe+f 

3 

As a result, 

n 

(l-\w\ 2 )-iG(w) = C^2 a ^ 1 -H 2 ) fc+f G fc+f H 



k=0 j 



II' I 



j fe=o 

= £fc>)(i-N a )~ 4 G(«i). 

3 

Now use estimate (|5.4|) . then Cauchy-Schwarz and finally ()5.2|) in the following way. 



2V-( 1 -I^| 2 ) 1 



J* |l-WjW| 

< c2(i-H 2 ) 2 |GK)| 2 E h _ ,, 

< c^a-H 2 )- 1 ^^)! 2 . 



Thus the factor (1 — |w| 2 ) 1 cancels and we have 

|GH| 2 <c2lGK)r 
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This last inequality gives at once 

J J \G(w)\ 2 (1 - \w\ 2 )- 2 dw < CJ2 \G{w 3 )Y 



or, in the upper half plane, 



Aj J \F(z)\ 2 S - 2 dxds <J2\F{z 3 )f 



In the wavelet notation this is 

r>oo poo 



poo poo 2 

A / W S S (i)f(-x,s) s~ 2 dxd S <Y / \Ws^)f(a J bm,a^) 

J J — oo 



Taking into account that S™(^) is admissible, we can apply (|3.ip and (|5.1[) to write 
the above inequality as 



(5.6) 



aii/ii 2 <e 



3,m 



f,Ta mJ D bm<j SZ(~) 



This is the left frame inequality. The sequence {T am -D\> m j S"(^)}j. m is thus a 
frame. □ 

In view of remark 5.1, the lattice result is a special case of our main theorem. 

Corollary 2. If b log a < then {T a j bm D ai 'S"(f)} iim is a frame of 

H 2 (U). 

Remark 5. Theorem 5.1 parallels theorem 3.1 in [12) . which states that, in 
the lattice case of the time- frequency plane, if the density of the lattice A is > n + 1 
(or if the size of A is < (n + l)^ 1 ) then the Gabor system {e 2mX2t H n {t — Ai) : 
A = (Ai,A2) £ A}, where H n stands for the Hermite function of order n, is a 
frame for L 2 (R). In particular, if one is dealing with the Von Neumann lattice 
with parameters a and b, the condition is ab < (n + 1) . Corollary 5.2 makes this 
analogy even more explicit. 

Remark 6. It is interesting to notice that, since Ber a f(z) = ^ t a e izt {Tf){t)dt, 
we clearly have 

k 



i k B( 



„fc+4 



/(*) = 



Ber 2 f(z). 



This allows to rewrite H5.5\) as 



^ (# >/(-*,.) = C£S^* + * (±) Ber^f(z) 



k=0 



which is, in some sense, reminiscent of Proposition 3.2 in |12j . 

Remark 7. Observe that combining the Paley-Wiener with the Plancherel the- 
orem, we have ||/||#2(u) = II/IIl 2 (o co) an< ^ ^ e resu ^ s °f Theorem 1 and 2 say also 
that we have a frame for L 2 (0, oo). 
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6. Further properties of the functions S"(t) 
The notation 



v 2 ^ (a)k(b)k_k 



F(a,b;c;x) = ^ 1 ^-x . 

k=o v Jk 

for the hypergeometric function is used in this section. Rewriting S" in this notation 
gives 

S%(t) = C a - n [ --it ) F(-n,- + l;a + l ;i ^ 



2 J v ' 2 ' ' i-it' 



(observe that the infinite sum becomes a polynomial of order n, since (— n)k = if 
k > 0). Composing the functions S"(t) with the fractional linear transformation 

2t - i 



2t 



(6.1) 
the result is 

(6-2) 5°(t)=r(| + l)(l-«)* +1 ff °(2r) 

where 

£( z ) = ^k F (_ n , " + 1; _ n _ Q /2 + 1; z) 

ni 2 

is a polynomial in z of degree n. This was pointed out in [29] . It was also shown 
that these polynomials satisfy the orthogonality 

dz 



(6.3) / g «( z )g%(z)\l-z\ a — =0 if m^n 

J\z\=l z 

and therefore are orthogonal on the unit circle with respect to the weight 

(6.4) w(z) = sin" -d6. 

This fact implies many properties, since there exists a very rich theory for orthogo- 
nal polynomials on the unit circle (see [30) and references therein and also chapter 
8 of |17j ). For example, the general theory assures that all the zeros of g"(z) lay 
within the unit disc. 

Remark 8. Setting a — — in Example 8.2.5 at [17] . and using the identity 
1 — e %e — 4 sin 2 | to write the measure (8.2.21) as {6.$ we recognize that the 
polynomials g"{z) are, up to a normalization, a family of orthogonal polynomials 
on the unit circle known as the circular Jacobi orthogonal polynomials. 

Remark 9. From \6. 3\) the functions z~^(l — z)^g^(z) are orthogonal on the 
circle and form a basis of the Hardy space on the unit disc 

H 2 (D) = {/ : / is analytic in D and sup / /(re 1 *) dt < oo} 

r<l Jo 

it is also clear that S% (t) are orthogonal on the real line ( the boundary of the upper 
half place). Since g"(z) are the circular Jacobi orthogonal polynomials, the basis 
functions z~^(l — z)% g"(z) are the circular Jacobi orthogonal functions and it is 
therefore natural to call S"(t) the rational Jacobi orthogonal functions. 
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Remark 10. From the general theory |17l (8.2.10)] follows that, if K n is the 
leading coefficient of the polynomial, then the sequence of polynomials {gn(z)} sat- 
isfies a three term recurrence relation 

K»ff£(0)ffn +l (*) + «*-lff£fl(<«-i(2) = Mn +l (0) + «n+lfl«(0)«]^(«) 

where K n is the leading coefficient of the polynomial. From the explicit representa- 
tion of the polynomials g%(z) it is easily seen that 

(f + l)n , a(f + l)n-l 
K n = j ; <Pn °) = 7T\ 

This three term recurrence relation provides a effective method for computational 
purposes: To evaluate the functions S"(t) it is sufficient to combine this recurrence 
relation with formulas &6.1\) and h6.Bjl . 
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